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ABSTRACT 
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Context. Gamma-ray burst emission in the prompt phase is often interpreted as synchrotron radiation from high-energy electrons 
accelerated in internal shocks. Fast synchrotron cooling of a power-law distribution of electrons leads to the prediction that the slope 
below the spectral peak has a photon index a = -3/2 {N(E) oc E a ). However, this differs significantly from the observed median 
value a as - 1. This discrepancy has been used to argue against this scenario. 

Aims. We quantify the influence of inverse Compton (IC) and adiabatic cooling on the low energy slope to understand whether these 
processes can reconcile the observed slopes with a synchrotron origin. 

Methods. We use a time-dependent code developed to calculate the GRB prompt emission within the internal shock model. The code 
follows both the shock dynamics and electron energy losses and can be used to generate lightcurves and spectra. We investigate the 
dependence of the low-energy slope on the parameters of the model and identify parameter regions that lead to values a > -3/2. 
Results. Values of a between -3/2 and -1 are reached when electrons suffer IC losses in the Klein-Nishina regime. This does not 
necessarily imply a strong IC component in the FermifLKI range (GeV) because scatterings are only moderately efficient. Steep 
slopes require that a large fraction (10-30%) of the dissipated energy is given to a small fraction (Sl%) of the electrons and that the 
magnetic field energy density fraction remains low (£ 0.1%). Values of a up to -2/3 can be obtained with relatively high radiative 
efficiencies (>50%) when adiabatic cooling is comparable with radiative cooling (marginally fast cooling). This requires collisions at 
small radii and/or with low magnetic fields. 

Conclusions. Amending the standard fast cooling scenario to account for IC cooling naturally leads to values a up to - 1 . Marginally 
fast cooling may also account for values of a up to -2/3, although the conditions required are more difficult to reach. About 20 % 
of GRBs show spectra with slopes a > -2/3. Other effects, not investigated here, such as a thermal component in the electron 
distribution or pair production by high-energy gamma-ray photons may further affect a. Still, the majority of observed GRB prompt 
spectra can be reconciled with a synchrotron origin, constraining the microphysics of mildly relativistic internal shocks. 

Key words, gamma-rays: bursts; shock-waves; radiation mechanisms: non-thermal 



1. Introduction 

The physical origin of the prompt emission in Gamma-Ray 
Bursts (hereafter GRBs) is still uncertain. The identification 
of the dominant energy reservoir in the relativistic out- 
flow, of the mechanism responsible for its extraction and 
of the processes by which the dissipated energy is eventu- 
ally radiated remains a major unresolved issue. There are 
three potential energy reservoirs : (i) thermal energy tha t 
can be radiated at the photo s phere (Meszaros & Rees I200C 



can be radiated at the photosphere (Meszaros & Keesl UUU 
Daig ne & Mochkovitchl l2002t iGiannios & Spruit! l2007t iPe^e 
2008; iBeloborodovl l2010t) . (ii) kinetic energy that can be 
extracted by shock waves propagating within the outflow 
and the n radiated by sh o ck-acc ele rated electrons (i nterna l 
shocks. iRees & Meszaros! dl994l) : iKobavashi et alJ dl997l) : 
Daigne & Mochko vitchl d 1998b ). or (iii) magnetic energy 
that can be diss i pated via the reconnec t ion of field lines 
dThompsonl 11994b iMeszaros & Reel 119971: ISoruit et all " 



Drenkhah n & Spruitl l200l iLvutikov & Blandfordl 



2001 



2003 



Giannios & Spruitl 120051) and then radiated by accelerated 
particles. In the two last cases, the expected dominant radiative 



Send offprint requests to: F. Daigne (daigne@iap . fr) 
* Institut Universitaire de France 



processes are synchrotron radiation and inverse Compton 
scattering. 

Observed GRB spectra can provide reliable constraints 
on the extraction mechanism and dominant radiative process. 
A typical GRB prompt emission spectru m is usually wel l 
described by a phenomenological model dBand et al.1 119931) 
where the photon flux follows N(E) oc E" at low energies and 
N(E) oc E 13 at high energies, with a smooth transition around 
£p, which is the peak energy of the vF v spectrum. Typical 

values in GRBs observed by BATSE are a 1 and fi 2.3 

for the low- and hig h-energy photon in dex and E p ~ 250 keV 
for the peak energy dPreece et alj |2000). It is well known that 
the observed value a ~ -1 is in clear contradiction with the 
predicted value for the synchrotron radiation from relativistic 
electrons (Pre ece et al.l 11998c iGhisellini et al.l 12000). One ex- 
pects a = -2/ 3 in slow coolin g regime and a = -3/2 in fast 
cooling regime dSarietall l998). the latter case being favored in 
the prompt GRB phase as the slow cooling regime would lead 
both to an energy crisi s and a difficulty to reproduce the shortest 
timescale variability dRees & Meszaros! 1 19941 : ISari et al1 [l996: 
IKobavashi etdl 19971) . 
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On the other hand, synchrotron radiation is a very natural 
expectation for the emission from shock-accelerated electrons. 
In GRBs especially, it is most probably at work in afterglows. 
Observations of prompt GRBs by the LAT instrument on 
board Fermi indicate that most GRBs do not show an addi- 
tional component at high energy (> lOO MeV) brighter or 
as bright as in the soft gamma-ray range (lAbdo et a 
lOmodei et alj|2009b . Prompt observations in the optical domain 
remain difficult but do not show strong evidence in favor 
of a bright additional component at lo w energy, with som e 
notable exceptions like GRB 0803 19B dRacusin et alj l2008h . 
This strongly favors synchrotron radiation compared to the 
synchrotron self-Compton (SSC) process for the emission 
observed in the soft gamma-ray domain (keV-MeV). Indeed, 
the latter requires some fine tuning to avoid a strong component 
either in the optical-UV-soft X-rays doma i n, or in the GeV 
range dBosnjak et alj 120091; IZou et all 120091: IPiran et al.ll2009b . 
Compared to SSC, synchrotron radiation has also a better 
ability - at least in the internal shock framework - to reproduce 
the observed spectral evolution : e.g. hardness-intensity and 
hardness-fluence correlations, evolution of the pulse width 
with energy channel, time lags (Daigne & Mochkovitch 1998; 
Ramirez-Ruiz & Fenimore 2000; Daigne & Mochkovitch 2003, 



Bosnjak 2010 in preparation). 

In this paper, we investigate a solution to steepen the low- 
energy slope of the synchrotron component and possibly recon- 
cile the synchrotron process with observed GRB prompt spec- 
tra. This solution is related to the steepening of the low-energy 
synchrotron slope by moderately efficient inverse Compton scat- 
terings i n Klein-Nishin a regim e, as suggested by_ Derishev et alJ 
(1200 lb: iBosniak et all d2009h : iNakar et alj (|2009j); see also 
iReesI d 1967b where this is mentionned in the general con- 
text of S SC radiation. It is an alternat i ve to the SSC scenario 
(see e.g. iPanaitescu & Meszarosl 120001 [Baring & Braby 2004; 
Kumar & McMahon 2008), to the comp tonization scenario 
dLiang et all 1 997t iGhisellini & Celottil 1999b . or to other propo- 
sitions to modify the standard synchrotron radiation, related 
to the timescale of the accele ration process dStern & Poutanen] 
l2004t lAsano & Terasawall2009b. the pitch-an gle distribution of 
electrons dLlovd-Ronning & Petrosianll2002b or the small scale 
struct ure of the magnetic field ( MedvedexfcOOOtlPe'er & Z hang 
2006). A summary of the measurements of the low-energy pho- 
ton index a and its distribution in GRBs is given in Sect. [2] Then 
Sect. [3] describes how the standard synchrotron spectrum is af- 
fected by additional processes such as inverse Compton scatter- 
ings or adiabatic cooling. It allows to identify physical condi- 
tions - in terms of intensity of the magnetic field, distribution of 
relativistic electrons, etc. - that lead to low-energy slopes steeper 
than the standard prediction a = -3/2. We discuss in Sect.|4]how 
such conditions could be found in GRB outflows. We compute 
expected pulse lightcurves and spectral evolution in the frame- 
work of the internal shock model and show that steep slopes are 
indeed expected in a large region of the parameter space. We 
summarize our conclusions and discuss future possible develop- 
ments in Sect. [5] 

2. A critical view on the observed distribution of 
GRB prompt spectral properties 

We examine here in detail the observed distribution of the low 
energy photon spectrum index a, as it provides a relevant crite- 
ria for the goodness of the emission model for GRB emission. 
To date the largest database of gamma-ray burst high time and 



energy resolution data was provided by the Burst and Transient 
Source Experiment (BATSE) (-20 k eV - 2 MeV) on boar d the 
Compton Gamma Ray Observatory. iKaneko et alj (|2006) pre- 
sented a systematic spectral analysis of 8459 time-resolved spec- 
tra from 350 GRBs (including 17 short events) observed by 
BATSE; this sample includes also gamma-ray bursts that were 
exa mined in previously publ ished catalogs of BATSE GRBs 
(e.g. lPreece et alj fl998, 2000). The reported distribution of a is 
apparently not c onsistent with the pred ictions of the simple syn- 
chrotron model: IKaneko et alJ |2006) showed that the median 
value for the time resolved spectra a = — 1.02^5'™ (long GRBs) 
and -0.87*q3q (short GRBs). The values obtained for the time 
integrated spectra indicate somewhat softer spectra, with the re- 
spective median indices —1.15^2™ and -0.99^2?^ for long and 
short events respectively. The slope a is distributed roughly sym- 
metrically around the median value. 

The results obtaine d by other ins t rumen ts are consistent with 
BATSE observations: Krimmetal. (2009) combined the Swift 
Burst Alert Telescope (BAT) and Suzaku Wide band All-Sky 
Monitor (WAM) data covering the broad energy band 15 to 
5000 ke V and report the distribution of a skewed toward slightly 
lower values, - 1 .23+0 .28 (for time-integrated spectra), while 
Pelangeon et al.l (|2008) for time integrated spectra of GRBs ob- 
served by High Energy Transient Explorer 2 ( HETE-2 ) in the en- 
ergy band 2-400 keV find a = -1.08±0.20. Similar results have 
been found recently by the Fermi Gamma-ray Space Telescope 
GBM and LAT in the broa d energy range ~ 8 keV t o > 100 GeV, 
e.g. the sample studied by iGhirlanda et all (1201 0b of 12 GRBs 
observed by Fermi displays various values for a, ranging from 
-1.26+0.04 (GRB 090618) to -0.55+0.07 (GRB 081222) for 
time-integrated spectra with determined peak energy. 

The observed distributions of GRB spectral parameters 
should however be considered with precaution. We point out the 
possible caveats for interpreting the results of spectral analysis: 

1. the spectral analysis are commonly performed on bright 
GRBs (with higher photon flux), which tend to have higher 
peak energies in general than dim GRBs ( M allozzi et al.1 
119951: iBor gonovo & Rvde 2001). For that reason the spec- 
tra with higher E pea k may be oversampled. The results of the 
time-resolved spectral analysis may be biased in the similar 
way: as the data were sampled more frequently during the in- 
tense episodes, the brighter portions of each burst may have 
more impact in the final distribution of spectral parameters 
( IKaneko et al.ll2006h : 

2. the low-energy photon spectra indices tend to correlate with 
the peak energy of the vF v spectrum, the sl ope a becom- 
ing softer when the peak energy is decreasing ( Kanek oet al.l 
2006; ICrider et al.ll 19971: iLlovd-Ronning & Petrosianll200a 
Ford et all 1 9951; iPreece et al.ll998b . This effect might be due 
to a combination of the curvature of the spectrum around the 
peak energy and the limited spectral energy range sampled 
by the instruments; 

3. as disc ussed bv IPreece et al.l ([l998) and iLlovd & Petr osian 
d2000b for the BATSE spectra, the data don't always ap- 
proach the GRB spectral low energy power law within the 
instrument energy range. If peak energy is close to the edge 
of the instrumental energy window, the low energy spectral 
power law may not have reached yet its asymptotic value. In 
that case lower va l ues of a are determined (i.e. softer spec- 
tra). |Ka£eko_eTaL| (|2006) attempted to account for this effect 
and applied as a better measure of the actual low energy be- 
havior the effective index a e /j for BATSE data, defined as 
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Fig. 1. Consequences on energetics of a steepening of the low-energy slope in the "Band" function. Left: the normalized "Band" 
function is plotted for B = -2.25 and a = -1.5 (standard value for synchrotron radiation in fast cooling regime) or a = -1 (mean 
value observed in GRB spectra). The dashed area corresponds to the fraction f e of the energy that should be removed by any process 
leading to such a change of the low-energy photon index. Here / e st 19%. Center: the fraction / e is now plotted as a function of B, 
assuming a low-energy photon index changing from a = - 1 .5 to 1 . Right: the fraction f e is now plotted as a function of a, when 
the photon index is changing from a = -1.5 to a. We assume B = -2.25. These three figures are plotted assuming that the process 
responsible for the change of the slope a does not affect the value of the break energy nor the high-energy tail of the spectrum. 



the tangential slope of the spectrum at 25 keV (the lower en- 
ergy limit of the BATSE window); 
4. the low-energy spectral index distributions of time- 
integrated and of time-resolved spectra are slightly different; 
it is expected due to the evolution of the spectral parame- 
ters during the integration time. The sharp spectral breaks 
are smeared over and the indices of time integrated spectra 
appear softer than i n the case of the time-resolved spectra 
dKaneko et alj |2006). These two last points could imply that 
the true low-energy spectral slopes in prompt GRBs are even 
steeper than the observed median value. 

Another important aspect to examine in the observed distri- 
butions of spectral indices concerns the contribution of an in- 
dividual GRB to the overall distribution. Since the brighter and 
longer events in general contribute with larger number of spectra 
in the BATSE spectroscopic catalog, we have computed distribu- 
tions of spectral parameters where each time bin of a given GRB 
is weighted by the corresponding fraction of the total fluence of 
the burst. In this way GRBs with different number of time re- 
solved spectra in their time histories hav e the same impact o n 
the overall distribution. Using the data bv lKaneko et alj (|2006), 
we find that: 

- Only 5% of GRBs have more than 50% of their spectra with 
very soft low energy slope, a < -1.5. Such slopes are prob- 
ably related to the spectral curvature around the peak energy 
and do not raise a problem for the standard synchrotron sce- 
nario; 

- 70% of GRBs in the sample have more than 50% of their 
spectra with a low energy photon index within the limits 
of the synchrotron model, -3/2 < a < -2/3 and 35% of 
GRBs have more than 50% of their spectra with a in the 
range -3/2 < a < -1; 

- Only 20% of GRBs in the sample have more than 50% of 
their spectra with a > -2/3. However most of these spectra 
have a close to -2/3. For instance, only 5% of GRBs in the 
sample have more than 50% of spectra with a > -0.3. 

From these values, it appears that the synchrotron radiation 
should in principle be compatible with most observed prompt 



GRB spectra. One should however understand why the mean 
value of a is steeper than the expected slope for the fast cool- 
ing regime, a = -3/2. 



3. Synchrotron radiation in radiatively efficient 
regime 

In this section, all quantities are given in the comoving frame of 
the emitting region. 

3.1. The standard prediction in fast cooling regime 

The synchrotron power of an electron with Lorentz factor y is 
given by 



(1) 



where B is the magnetic field. If the source is relativistically 
expanding, as expected in gamma-ray bursts, adiabatic cooling 
competes with radiation. This cooling process occurs on a typ- 
ical dynamical timescale f ex . It is then convenient to define F c 
as the Lorentz facto r of electrons w hose synchrotron radiative 
timescale equals f ex (Sari et al. 1998J) : 



r, = 



6nm e c 



(Tj 



B 2 U 



(2) 



If the GRB prompt emission comes from the radiation of rel- 
ativistic electrons, it is necessary, both to allow for the short- 
est timescales observed in GRB lightcurves and to minimize the 
constraint on the energy budget, that these electrons are radia- 
tively _e^oi«U 2 J^e 1 _flw^heir radiative t i mesca l e is shorter than 
U (jRees & Meszarosl [19941 ISari et alJ [19961 iKobavashi et all 
1 19971) 7 When only synchrotron radiation is considered, this is 
equivalent to the condition that most injected electrons by the ac- 
celeration process must have y > F c . The resulti ng synchrotron 
spectrum has been described bv lSari et alj d 1998b when the ini- 
tial distribution of relativistic electrons is a power-law of slope 



4 



Daigne, Bosnjak & Dubus: Reconciling observed GRB prompt spectra with synchrotron radiation ? 





Fig. 2. The effect of inverse Compton scatterings in Klein-Nishina regime on the fast cooling synchrotron spectrum. The 

normalized synchrotron spectrum defined by v w v | syn lrif c T m m e c 2 is plotted as a function of the normalized frequency v/v m , as well 
as the corresponding photon index dlnu v /dlnv-l. All spectra in thick solid line are computed numerically using a detailed radiative 
code including synchrotron radiation and inverse Compton scatterings (see text). Spectra in thin dotted line are computed with 
synchrotron radiation only. For clarity purposes, all other processes (adiabatic cooling, synchrotron self-absorption, yy annihilation) 
are neglected. A ratio r c /r m = 10~ 6 is assumed to ensure that all electrons are in fast cooling regime. The maximum Lorentz factor 
of electrons is fixed to F max = 10 4 r m ; n so the high energy cutoff in the synchrotron spectrum appears at the same frequency in 
all cases plotted here. The four panels correspond to increasing values of the w m parameter, w m = 0.01, 1, 100 and 10 4 from 
the top-left to the bottom-right panel, i.e. to a growing importance of Klein-Nishina corrections for inverse Compton scatterings. 
In each panel, the six curves in solid line correspond to increasing values of the Y parameter, Y = 0.1, 1, 10, 10 2 , 10 3 and 10 4 , 
i.e. to a growing efficiency of inverse Compton scatterings. The table inserted in each panel lists the values of the ratio £i c /£ syn 
of the inverse Compton component (not plotted here) over the synchrotron component. In the bottom-left panel (w m = 100), the 
synchrotron spectrum obtained assuming a slow injection over f ex (see text) is also plotted in dashed line for comparison. 



—p above a minimum Lorentz factor F m . If synchrotron self- 
absorption is neglected, three asymptotic branches are predicted 

(£) if v c <v<v m , (3) 



'Isyn 



nf c r m m e c 2 



where u v is the final photon energy density at frequency v and 
nf c is the initial density of relativistic electrons. The break fre- 
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quency v m (resp. v c ) is defined as the synchrotron frequency for 
an electron with Lorentz factor T m (resp. T c ). This asymptotic 
spectrum shows clearly that the predicted photon spectral slope 
a below the peak of vF v is -3/2, in apparent contradition with 
observations. 



3.2. The effect of Inverse Compton scatterings in 
Klein-Nishina regime 

The mean value of a observed in the BATSE spectroscopic cat- 
alog (see Sect.[2ji is close to -1. As seen in Fig.[T| for typical 
values of the high-energy photon index /3 between -2 and -3, 
a change of the low-energy slope of the Band function from 
a = -1.5 to -1 requires only a sub-dominant radiative pro- 
cess that can transfer about 20 - 40 % of the energy from the 
synchrotron component into another component. A natural can- 
didate is inverse Compton scattering of synchrotron photons by 
relativistic electrons. Indeed, this process necessarily takes place 
in the emitting region. Two parameters can be introduced to eval- 
uate the importance of inverse Compton scatterings. The first pa- 
rameter, w m , measures if scatterings occur mostly in Thomson 
regime (w m <K 1) or if Klein-Nishina corrections are important. 
It is defined by 



Tin 6 n 



(4) 



The second parameter, Fxh, measures the intensity of the inverse 
Compton process in the Thomson regime and is defined by 



Tni 



cr T nf c ct. 



syn (T m ) X ^ T^ 



where 



hyn (7) 



ym e c 

^syn(r) 
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(5) 



(6) 



is the synchrotron timescale of an electron with Lorentz factor y. 
In Eq.|5]the first term, cr x «g CC cf syn , is the Thomson optical depth 
associated with relativistic electrons in fast cooling regime, 
and the second term, |F^ , corresponds to the typical boost of 
a photon scattered by an electron at Lorentz factor F m , if the 
scattering occurs in Thomson regime (w m <K 1). Therefore, if 
w m <K 1 and ixh ^ 1 . the ratio of the total energy & K radiated 
in the inverse Compton component over the total energy £ syn 
radiated in the synchrotron component is simply given by 
£j c /£ syn Fxh- Still in the Thomson regime (yv m <K 1), if Ixh 
is large, the effective radiative timescale becomes ~ f syn /Fxh 
and fi; c /fi syn ^ y/Yjh- Finally, when Klein-Nishina corrections 
are important (w m <: 1), both the cross section and the boost in 
frequency are signif icantly reduced so that £j c /£ syn <sc Yxh (see 
Bosnjak et al. 2009, for details). 

If the magnetic energy density represents a fraction of the 
local energy density in the emitting region, and if the energy in- 
jected into accelerated relativistic electrons represents a fraction 
e e of the same energy reservoir, Eq.[5]can be simplified to give 



Fni - 



p-2 e e 
p-l €b 



(7) 



When considering only synchrotron radiation and inverse 
Compton scatterings, the spectral shape of the synchrotron com- 
ponent, i.e. vM v | syn /nl cc F m m e c 2 as a function of v/v m , depends 
only on these two parameters w m and Fph, in the limit of ex- 
treme fast cooling (F c «: r m ). It is well known that in Thomson 



regime (w m <K 1), as the electron cooling rate due to inverse 
Compton scatterings remains proportional to y 2 like for the syn- 
chrotron power, the spectral shape of the sy nchrotron compo- 
nent is un-affected b y the scatterings (see e.g. lSari & Esin 2001; 
Bos njak et alj2009h . To change the value of a, the physical con- 
ditions in the emitting regions must necessarily be such that 
w m ^ 1. Then, Klein-Nishina corrections become important 
for most of the scatterings and the dependence of the el ectron 
cooling rate on y differs from y 2 . IDerishev et alj ((2001) have 
shown that it results in a steeper slope a, that can potentially 
reach the value a - -I. We have shown in a previous paper that 
this change of the synchrotron slope due to inverse Compton 
scatterings in Klein-Nishina r egime was indeed observed in de- 
tailed radiat i ve calc ulations dBosniak et alj2 009). More recently, 
iNakar etal] ([2009) have presented a complete analytical esti- 
mate of the asymptotic synchrotron spectrum in the presence of 
inverse Compton scatterings, and have shown that the asymp- 
totic slope a = -1 is expected in the synchrotron component 
below the peak of vF v wherQ F c <K F m (fast cooling), w m » 1 
(Klein-Nishina regime) and 

< 3 < Y < w 3 m . (8) 

Therefore t h e ana lytical work of IDerishev et al.1 (l200ll) and 
INakar etal] ((2009) and the numerical study of Bosnja k et al.l 
(2009) indicate that the effect of inverse Compton scatterings in 
Klein-Nishina regime on the synchrotron component seems to 
be a promising possibility to explain observed values of a in the 
range [-3/2; -1]. However, as all possible values in this interval 
can be observed, it is necessary to go a step further than analyti- 
cal estimates based on asymptotic behaviours. In Fig. [2] we show 
the evolution of the synchrotron spectrum when Fxh is increas- 
ing, for different values of w m (and for F c <s F m ). The spectra 
are co mputed using the radiative code described in Bosnjak et al. 
(2009) that solves simultaneously the equations of the time- 
evolution of the electron and photon distributions and that in- 
cludes most relevant processes (adiabatic cooling, synchrotron 
radiation and self-absorption, inverse Compton scatterings and 
photon-photon annihilation). To focus on the effect described 
in this subsection, the spectra in Fig. [2] are computed including 
only synchrotron radiation and inverse Compton scatterings and 
do not take into account the other processes, whose impact is dis- 
cussed later in the paper. The slope a is clearly steepening con- 
tinuously from -3/2 to -1. In Fig. [3] we have color-coded in the 
diagram w m — iTh the value of the low-energy photon index of the 
synchrotron spectrum. This slope is only a representative value 
as the synchrotron spectrum below the peak shows some curva- 
ture with an evolving slope (see Fig. [2]). In practice, we adopt 
for a the clear maximum of the slope below v m . If a — -3/2 
appears indeed universal when inverse Compton scatterings are 
negligible (Fxh ^ 1) or Klein-Nishina corrections unimportant 
(w m <sc 1), the situation is rather different in the quarter of the di- 
agram where w m > 1 and Fxh > 1 . In particular, a region where a 
evolves from -3/2 to - 1 is well defined, whose boundaries agree 
reasonably well with the condition given by Eq.[8] In the same 
figure, lines of constant ratio £i c /£ S yn are also plotted. In agree- 
ment with the analysis made at the beginning of this subsection, 
this ratio is typically in the range 0.1-1 in the region of inter- 
est (it is of course much smaller than Fph due to Klein-Nishina 



1 In lNakar et alj i2009h . the situation where the slope a = - 1 is pos- 
sible corresponds to cases lib and lie, which - in addition to fast cooling 

- are limited by ir m /r m ) < e c / f B < (r ra /r ra 1 , where the authors de- 
fine y = m c c 2 //iy syn (y) so that w m = r m /r m , leading to the conditions 
given in Eq.[8] 
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Fig. 3. The low-energy slope a of the fast cooling synchrotron 
spectrum in the presence of inverse Compton scatterings in 
Klein-Nishina regime. In the w m -Y T h plane, the value of the 
photon index a of the synchrotron spectrum below the peak of 
vF v is color-coded. All spectra have been computed numerically 
using a detailed radiative code including synchrotron radiation 
and inverse Compton scatterings (see text). Black thick lines of 
constant ratio Li c /L syn are plotted on top of this diagram. Cases 
where L; c /L syn > 5 in the top-left corner are not plotted as we 
focus in this paper on situations where the synchrotron compo- 
nent is dominant. In addition, thin dashed l ines s how the location 
of regions lib and lie from Na kar et alJ J2009), where a slope 
a — > - 1 is predicted. 



corrections to the cross section and the boost in frequency). It is 
important to note that a - 1 can be reached in practice for not 
too extreme values of the parameters, typically w m = 100 - 10 4 
and Y ^ w m . 

The maximum value of a that can be reached depends 
on the assumptions about the electron acceleration process. 
In principle, one would wish to follow the full 'magneto- 
hydrodynamical' evolution of the shocked region at the plasma 
scale. The framework used in the present study to follow the 
dynamics of shocks in a relativistic outflow limits us to make 
simple assumptions on the electron injection timescale. Two ex- 
treme cases are possible. If electrons are injected reg ularly over 



a tim escale f; n j ec comparable to f ex (as assumed in iNakar et aH 



2009), the slope will never be steeper than a - -I and in prac- 
tice will be a little less steep than this absolute limit for reason- 
able choice of parame ters (not too extreme w m and Y parameters, 
see lNakar et alJl2009l and bottom-left panel in Fig.|2). If the in- 
jection occurs faster, the limit a = -1 is more easily reached 
and even steeper values can be obtained for extreme sets of pa- 
rameters (see the maximum slopes obtained for F t h = 10 4 and 
w m = 100 or 10000 in Fig.[2j». Calculations presented in this pa- 
per corresponds to the regime where f; n j ec <s t ex . 



3.3. Additional effects 
3.3.1. Adiabatic cooling 

The process described in the previous subsection offers a physi- 
cal interpretation of the observed values of the low-energy pho- 
ton index in the range [-3/2; -1]. However the value a = — 1 
does not appear as a limit in the observed distribution of a and 
an additional explanation has to be found for the steeper slopes. 
When electrons are in slow cooling regime (T c » r m ), the 
predicted value of the p hoton index below the peak of vF v is 
a - -2/3 dSari et al.lll998l) . Therefore, it is now tempting to as- 
sociate spectra with -1 < a < -2/3 to this regime. There is 
however a problem due to the radiative efficiency / ra d = u y luf c , 
where uf c is the initial energy density injected in relativistic 
electrons and u y the final energy density of the radiated photons. 
In the slow cooling regime / ra d is low, which increases the re- 
quired energy budget to an uncomfortable level. Here, we rather 
consider the situation where electrons are in fast cooling regime 
but not deeply in this regime, i.e. F c r m rather than T c <sz F m 
("marginally fast cooling regime"). 

To illustrate this situation we plot in Fig. |4] the evolution of 
the synchrotron spectrum for an increasing ratio r c /r m and for 
different values of (w m ; Yjh) representative of the different re- 
gions in the diagram of Fig. [3] As expected a break in the spec- 
trum appears at frequency v Ci eff , i.e. at the synchrotron frequency 
of electrons with Lorentz factor r c e jf, whose radiative timescale 
is equal to the dynamical timescale f ex . We have r c>e ff < F c due 
to inverse Compton scatterings. The photon index below v c?e ff is 
-2/3. Therefore, when v c>e s is close to v m the observed photon 
index can be very close to this asymptotic value, even in fast 
cooling regime. This is well seen in Fig. [4] for w m = 100 and 
y Th = 100 (top-right panel) or w m = 10 4 and Y Th = 10 4 (bottom- 
right panel) and for r c /F m = 1 (/ lac j 0.6 - 0.7 in this case). 
When the efficiency of inverse Compton scatterings is reduced, 
v c , e ff is closer to v c and the same effect can be seen for lower val- 
ues of the ratio F c /r m . This is the case for instance for w m - 0.01 
and Irh = 0.1 (top-left panel) or for w m = 10 4 and irh = 100 
(bottom-left panel) and for r c /r m = 0.1-1 (/ rad =* 0.6-0.7 -> 1 
in this case). 

We have plotted in Fig.|5]the same diagram as in Fig. [3] now 
including the effect of adiabatic cooling for different values of 
the ratio r c /r m . The representative value of a is selected in 
the same way as in Fig. [3] but the maximum is not always as 
clearly defined as in the r c /r m = case (see for instance the 
curves for r c /F m = 0.1 in the four panels of Fig.|4]i. Clearly, 
when the intermediate maximum of the spectral slope below 
v m disappears and a is given the asymptotic value -2/3 (green 
region in Fig. [5), the comparison with observations of the 
low-energy photon index becomes difficult : depending on the 
location of the peak energy E p , and of the low-energy threshold 
of the instrument, any value of a between -2 and -2/3 can be 
measured. 

The diagrams in Fig. [5] show that already for r c /r m = 0.1, 
the asymptotic slope a = -2/3 is reached in a large region of 
the w m - Fxh plane, together with a large radiative efficiency. 
Even for F c /r m = 1, a = -2/3 is observed in a large region 
where the radiative efficiency is still larger than 66%. It is only 
for r c /r m = 10 that the slow cooling regime dominates the 
diagram, high radiative efficiency being found together with 
a — -2/3 in only a very small area. 
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Fig. 4. The effect of adiabatic cooling on the fast cooling synchrotron spectrum in presence of inverse Compton scatterings. 

The normalized synchrotron spectrum defined by v M v | syn lrif c Y m m e c 2 is plotted as a function of the normalized frequency v/v m , 
as well as the corresponding photon index d In u v /d In v - 1. All spectra in thick solid line are computed numerically using a 
detailed radiative code including synchrotron radiation, inverse Compton scatterings and adiabatic cooling (see text). Spectra in thin 
dotted line are computed without inverse Compton scatterings. All other processes (synchrotron self-absorption, yy annihilation) are 
neglected. Each panel corresponds to a different set of parameters (w m , irh) indicated in the top-left corner. In each panel, spectra 
are plotted for increasing ratios r c /r m = 0.01, 0.1 , 1 and 10. As in Fig. [2] the maximum Lorentz factor of electrons is computed 
with a fixed ratio r max /r m i n = 10 4 . Spectra in dashed lines are radiatively inefficient (slow cooling regime). The table inserted in 
each panel lists the values of the radiative efficiency of the electrons. 



3.3.2. Other effects 

Other processes may influence the final spectral shape. 
Photon-photon annihilation produces a cutoff at high energy 
dGranot et al.l l2008: Bos niak et al.l l2009). For all cases presented 
in this paper, we have checked that the opacity for this process 
was extremely low below ~ 100 Me V. Photon photon annihi- 



lation can affect the tail of the inverse Compton component 
at high-energies. The fraction of the radiated energy which is 
reinjected in pairs via yy — > <? + e~ is usually, but not always, 
small in the examples shown in Sect. [4] For example, in the 
three cases defined in Sect.[4] it is typically less than 10~ 3 
in case C, less than 0.05 in case A and between 0.1 and 0.3 
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Fig. 5. The low-energy slope a of the synchrotron spectrum in the presence of inverse Compton scatterings in Klein-Nishina 
regime, including the effect of adiabatic cooling. Same as Fig. [3] now including adiabatic cooling for r c /T m = 10~ 2 , 10 _I , 1 
and 10. Black solid lines of constant radiative efficiency f m & are plotted on top of the three last diagrams. In the first panel, the 
radiative efficiency is always close to 100 %. For larger ratio F c /r m , the radiative efficiency remains larger than respectively 96 % 
(F c /F m = 10- 1 ), 65 % (F c /F m = 1) and 28 % (T c /r m = 10). 



in case B. For collisions occuring at low radius or with low 
Lorentz factors, yy annihilation could be even more important. 
When the fraction of the energy in annihilated photons is non 
negligible, the resulting radiation of the created pairs could 
affect the spectral shape even at low energy, and modify the 
results presented here. Despite its potentially interesting impact, 
we defer to future work the investigation of such cases, due to 
the additional complexity it involves for the computation of the 
radiated spectrum. Numerical approaches to solve such a highly 
non linear problem including thermalization effects have been 



proposed by iPe'er & Waxmanl d2005l); lAsano & Inoud d2007l) : 
iBelmont et all (12008b : IVurm & Poutanenl (120091) ~ 



At low energy, the synchrotron self-absorption can also 
steepen the spectrum. This effect is included in simulations pre- 
sented in Sect. [4] and is always negligible in the soft gamma- 
ray domain, in agreement with the standard predictions for 
the synchrotron fast cooling regime. Indeed the timescale for 
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self-absorption at y r 
iBosnjak et ai1l2009l) 



is given by (see e.g. equation (28) in 



fa (V m ) 4tTV^ 



UkeV/ lis/ 



10 



-6 



(9) 



where 1 keV is taken for a typical value of the peak energy 
in the comoving frame and other parameters are given rep- 
resentative values for internal shocks. The timescale for self- 
absorption around v m is therefore always much larger than all 
other timescales (dynamical or radiative) and the self-absorption 
process is negligible in the soft gamma-ray range. 

In addition to the details of the radiative processes, the pre- 
cise shape of the electron distribution can also have an impact 
on the final spectrum. Here, we assume a power-law distribu- 
tion. More complex distributions showing several components 
(e.g. Maxwellian distribution + non-thermal tail) are observed 
in some s imulations of parti c le acceleration in u ltra-relativistic 
shocks (ISpitkovskvll2008blal: iMartins et al1l2009l) . Such results 
would need to be confirmed for the mildly relativistic regime 
of interest for the prompt GRB emissio n. In the ultra-relativistic 
regime relevant for the afterglow, Giannios & Spitkovskvl fcOOg) 
have shown that the Maxw ellian component could have an ob- 
servable signature. However lBaring & Brabvl d2004l) find that the 
non-thermal electron population should dominate in the prompt 
phase. We leave to a future work the study of the consequences 
of more complex electron distributions on the observed GRB 
prompt spectra. 



4. Constraints on the internal shock model 

4.1. General constraints on the physical conditions in the 
emitting regions 

We have shown in Sect. [3] that the spectrum resulting from 
synchrotron radiation in the presence of inverse Compton 
scatterings in Klein-Nishina regime can account for observed 
low-energy photon index a = -3/2 to -1, and that the additional 
effect of adiabatic cooling with T c 55 F m can lead to steeper 
slopes up to a = -2/3 in a regime where the radiative efficiency 
is still reasonably high (f m d ^ 50 %). In principle, this allows 
to reconcile the synchrotron process with the observed spectral 
parameters in most GRB spectra (see Sect.|2}. However it is still 
necessary to demonstrate that the physical conditions identified 
in Sect.[3]can be reached in the emitting regions within GRB 
outflows. These conditions are approximatively w m <: 1 and 
Wm 3 < T-rh ~ w 3 m for changing the slope from a = -3/2 to —1 
with the effect of IC scatterings in Klein-Nishina regime, and 
r m /r c — 0.1 - 1 to get steeper slopes up to a = -2/3 with the 
effect of adiabatic cooling. 



We assume that the prompt gamma-ray emission is produced 
in a relativistic outflow ejected by a source at redshift z. We con- 
sider an emitting region at radius R within the outflow, with a 
Lorentz factor F». We do not specify at this stage the physical 
mechanism responsible for the energy dissipation in this region, 
leading to the presence of a magnetic field B and a population of 
relativistic electrons with a minimum Lorentz factor F m . We as- 
sume that the medium is optically thin, i.e. tt = o-jnl cc ct ex < 1. 
The observed peak of the synchrotron spectrum is given by 

(1 +z)hv mobs * 380keV (— \(^\ (— — ) . (10) 
v ' ' \300/\10 4 / V3000G/ 
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Fig. 6. Exploration of the internal shock parameter space (1). 

All cases fulfilling the conditions of (i) transparency; (ii) radia- 
tive efficiency; and (iii) synchrotron dominance at low-energy 
(see text) are plotted with black dots in three planes : £p, b s 
vs a (top-left panel); £; c /£syn vs a (bottom-left panel); £p, b s 
vs £j c /£ syn (top-right panel). The ratio £; c /£ syn does not take 
into account the fraction of high energy photons that are sup- 
pressed by yy annihilation. In addition, the cases where the syn- 
chrotron spectrum peaks in the gamma-ray range (lOkeV < 
£p,ot>s i= lOMeV) and the low-energy photon index is in the 
range -3/2 < a < — 1 (resp. -1 < a < -2/3) are plotted in 
red (resp. blue). 



The dynamical timescale relevant for adiabatic cooling can be 
estimated by 



-l 



,„ =Lls ( * ) -L. 

\10 13 cm/\300 

From Eq.|2] this leads to 



B 



.3000G; 



) 10 4 



(ID 



(12) 



The two parameters governing inverse Compton scatterings are 
given by Eq.|4]and Eq.|5]: 



ioo(H(^) 

\10 4 / V3000G/ 



(13) 



and 



Lorentz factors F* above 100 are necessary to avoid the pres- 
ence of a high-ene rgy cutoff in the spect rum due to yy an- 
nihimation (see e.g. Lithwick & Sarill200ll) . Higher values are 
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Fig. 7. Exploration of the internal shock parameter space (2). The distribution of the three parameters defining the relative 
importance of inverse Compton scatterings, Klein-Nishina corrections and adiabatic cooling, F t h, w m and r c /T m are plotted for 
all cases presented in Fig. [6] In addition the distributions of the same parameters for cases where the synchrotron spectrum peaks 
in the gamma-ray range (lOkeV < E po \, s < lOMeV) and the low-energy photon index is in the correct range are also plotted 
(-3/2 < a < -1: red; -1 < or < -2/3: blue). 



even required in some bursts detected by Fermi/hA^T. For in- 
stance F ^ 600 - 900 has bee n derived by the F ermilLKY col- 
laboration for GRB 080916C dAbdo et alj 12009b . Radii in the 
range 10 13 -10 15 cm are expected as the prompt GRB emission 
should mainly occur above the photospheric radius and below 
the deceleration radius. Then, the main constraint comes from 
the fact that in the proposed scenario gamma-ray photons must 
be produced directly by synchrotron emission. From Eq.[l0l this 
is always possible if electrons can be accelerated to very high 
Lorentz factors. Then Eq.Qj] and (TBI) show that the physical 
conditions listed above and leading to -3/2 < a < -1 can in- 
deed be reached in GRBs. In addition Eq.[T2l indicates that the 
"marginally fast cooling regime" leading to a — > -2/3 can in 
principle also be reached, especially for emission at larger ra- 
dius, where the magnetic field could be lower. 

4.2. Impact on the microphysics parameters in internal 
shocks 

The internal shock model allows a self-consistent evaluation of 
the physical conditions (i.e. T*, R, B, T m , etc.) for each emitting 
region. It is then possible to identify the pertinent range of the 
parameters leading to steep low-energy slopes : the properties 
of the relativistic outflow (Lorentz factor, kinetic energy) and 
the parameters describing the microphysics at work in shocked 
regions (particle acceleration, magnetic field amplification). 

We consider first collisions between t wo equal-mass shells 
dBarraud et al.1 l2005t iBosniak etaD 120091) . More realistic out- 
flows are considered in the next subsection. The parameters 
defining the dynamical properties of a collision are the mean 
Lorentz factor in the outflow f , the ratio k of the Lorentz factor 
of the faster shell over the Lorentz factor of the slower shell, 
the time separation t between the ejection of the two shells 
and the kinetic energy flux E injected in the outflow. These 4 
parameters allow to estimate the radius of the collision, as well 
as the physical conditions in the shocked region (Lorentz factor 
r„, comoving density p„ and comoving specific energy density 
e*). In addition, four microphysics parameters are necessary 
to estimate the distribution of relativistic electrons and the 
magnetic field : and e e are the fraction of the energy density 



p,e, that is injected in the magnetic field and the relativistic 
electrons, respectively. We assume that only a fraction ( of the 
available electrons are accelerated in a non-thermal distribution 
and that this distribution is a power-law with index —p. The fu ll 
description of this model can be found in Bosni ak~et al.l J2009). 

We explore the parameter space of this model, assuming 
a constant value e e = 1/3 and a constant electron slope 
p = 2.5. A high value of e e seems unavoidable in internal 
shocks to maintain a reasonable efficiency of the process. We 
have checked that our results are not affected much by taking 
e e = 0.1. As we are mainly interested in the low-energy photon 
index a, the choice of p is not very important. The present value 
p = 2.5, leads to a high-energy photon index B -2.25 close 
to the mean value observed in BATSE spectroscopic catalog 
dPreece et al.l 120 001 We then compute the observed spectrum 
of 50400 collisions for logf =1.5, 2, 2.5 and 3; k =2.5, 5, 
7.5 and 10; logr =-2, -1, 0, 1 and 2; log£ =50, 51, 52, 53, 
54 and 55; log£ =-4, -3, -2, -1 and 0; loge B = -5.5 -> -0.5 
every 0.25. These spectra are c omp uted with the radiative 
code developed by Bos niak et alJ d2009h . including all relevant 
processes : synchrotron radiation, inverse Compton scatterings 
and adiabatic cooling, as considered in the previous section, but 
also synchrotron self-absorption and yy annihilation. 

We keep only cases which fulfill the following conditions : 
(i) the shocked region is transparent (cr T « ± c? e x < 0.1, where 
« ± is the final lepton density in the shocked region, taking into 
account pairs that wer e created by yy annih ilation but neglecting 
pair annihilation, see Bosni aket alJ 12009) ; (ii) electrons are 
radiatively efficient (/ rac ) > 0.5) ; (iii) synchrotron radiation 
is dominant at low-energy (M v | syn > 10 u v \ ic at the frequency 
of the synchrotron peak). This last condition eliminates a few 
cases where inverse Compton scatterings are so efficient that the 
synchrotron component is hardly observed. After this selection, 
~ 75% of cases are suppressed, because of a too large optical 
depth (condition (i): ~ 50% of cases), a too low radiative 
efficiency (condition (ii): ~ 20% of cases) or a negligible 
synchrotron emission (condition (iii): ~ 10% of cases). For 
each spectrum, assuming a source redshift z = 1, we compute 
the observed peak energy of the synchrotron component £p j0 bs, 
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the low-energy photon index a below the peak, and the ratio 
£ic/£.syn of the inverse Compton component over the syn- 
chrotron component. All models fulfilling the three conditions 
listed above are plotted in Fig. [6] This figure illustrates that the 
internal shock model with a dominant synchrotron process in 
the soft gamma-ray range (BATSE, Fermi/GSM) allows a large 
range oflow-energy photon index between -3/2 and -2/3 anda 
large range of peak energies (including very high peak energies 
as observe d in s ome bright bursts such as G RB 0809 16C, 
lAbdo et al.1 12009L for which IWang et all 120091 have recently 
shown that the slope a -1 can be explained by the effect of 
IC scatterings in KN regime on the synchrotron spectrum; and 
very low peak energies as expected in X-ray rich gamma-ray 
bursts or X-ray flashes, Heis e et al. l l200ll:ISakamoto et al]l2005l 
2008). All these cases have / ra d > 0.5 and even _f ra d > 0.9 in half 
of the cases. For cases with lOkeV < iip^bs ^ lOMeV, the ratio 
£ic/£syn is typically in the range 10~ 3 -1. This is in agreement 
with the indication from the Fermi/LAI GRB detection rate 
that most GRBs do not have a st rong additional c omponent 
between 100 MeV and 10 GeV dGranot et all 1201 Oal) . Note 
that the density of points in Fig. [6] has no physical meaning as 
the distribution of the physical parameters in GRB outflows 
is unknown. This figure only illustrates the range of observed 
values that can be expected in the internal shock model. 

We plot in Fig.|7]the distributions of w m , Fxh and r c /r m for 
the same models. The values of the low-energy photon index a 
are in full agreement with the analysis made in section Sect. [3] 
How are such values obtained ? The distributions of dynamical 
(T, k, t and E) and microphysics and es) parameters for all 
models that fulfill the three conditions listed above and have in 
addition a synchrotron spectrum that peaks in gamma-rays with 
a low-energy index in the expected range (-3/2 < a < -2/3) 
show that low-values of ( are needed to produce gamma-rays 
(typically £ = 10~ 3 - 10~ 2 ), and that low values of 6b favor 
steeper low-energy photon indexes. 



4.3. A few examples of synthetic GRBs 

In a more realistic description of the internal shock model, each 
pulse in GRBs with complex multi-pulses lightcurves is associ- 
ated with the propagation of a shock wav e within the relativistic 
outflow (IDaigne & Mochkovitchl 120001: iMimica et all l2007t 
iMimica & Alovl 120101) . This propagation implies an evolution 
of the physical quantities in the shocked region, especially the 
density and therefore the magnetic field. This leads to a spectral 
evolution within each pulse that h a s already be e n partially de- 
scribed in IDaigne & Mochkovitchl (Il998l 120031) : iBosniak et al.1 
(120091) . The model developed bv IBosniak et ail (120091) couples 
a dynamical simulation of internal shock propagation within a 
relativistic outflow, and a detailed radiative code. This allows 
to predict the lightcurves and spectral evolution in pulses for 
different assumptions regarding the physical conditions in 
the outflow. In order to compare the results with observed 
distributions of spectral parameters, we face difficulties due 
to several possible biases, as described in Sect. [2] To make a 
full comparison, one should generate noise in our synthetic 
bursts and take into account the response function of a given 
instrument before fitting the resulting spectrum by a Band 
function. We did not follow this procedure as our primary goal 
is to identify the theoretical limits for the prediction of the 
low-energy slope. We computed theoretical spectra over time 
bins of duration 0.25 s and measure the slope below the peak 
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Fig. 8. Case A : an example of a pulse generated by an in- 
ternal shock with synchrotron radiation in pure fast cool- 
ing. The top panel plots the evolution of the three main pa- 
rameters shaping the radiated spectrum (see Sect.[3} : w m , Tni 
and F c /F m . The four other panels show the lightcurves in dif- 
ferent energy channels corresponding to the GBM and the LAT 
on board Fermi. The respective contributions of synchrotron and 
inverse Compton are also indicated. In the two GBM panels, the 
evolution of the low-energy slope and peak energy of the soft 
gamma-ray component is also plotted. In this case, the standard 
fast cooling synchrotron slope a = -3/2 is found for the whole 
duration. 



energy in the same way as in Sect.[3] 

In the examp l es belo w, we adopt the same reference case as 
in Bosni ak et al.1 d2009l) : a single pulse burst is generated by 
a relativistic ejection lasting for f w = 2 s with a constant E 
an d a Lorentz factor in creasing from 100 to 400 (see Figure 1 
in IBosniak et aD l2009h . Constant microphysics parameters are 
assumed during the whole evolution. This is a simplifying as- 
sumption due to our poor knowledge of the physical processes 
at work in mildly relativistic shocks. However, as the diversity of 
GRBs and their afterglows seem to indicate that these parameters 
are not universal, they are most probably evolving with shock 
condit ions, which could impact the spectral evolution within a 
pulse (Daigne & Mochkovitch 2003). We adopt here e e = 1/3 
and p = 2.5 and adjust f and £b to have the peak energy of the 
pulse well within the GBM range. We consider the following 
examples to illustrate the possible range of a : 

- CaseA:£'= 10 54 erg/s, e B = 1/3 and ( = 3xl0~ 3 . This case 
is plotted in Fig. [8] and shows the standard slope a = -3/2. 
The peak energy is £ P;0 bs — 800 ke V at the peak of the pulse. 

- Case B : E = 10 54 erg/s, e B = 10~ 3 and £ = 10~ 3 . This case 
is plotted in Fig. |9] and shows a steeper or =a — 1, The peak 
energy is £ p , bs — 700 keV at the peak of the pulse. 

Intermediate values of 6b between case A and B would lead 
to intermediate values of a between -3/2 and -1. In both 
examples, it appears clearly that all spectral parameters are 
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evolving during a given pulse. The evolution for the peak energy 
is stronger in case A than in case B, whereas the low-energy 
photon index a evolves more strongly in the second case. Note 
that in these examples, the direct emission from pulse ends at 
fobs ~ 5.5 s. The flux observed after this time is due to the high 
latitude emission. In more complex bursts with multi-pulses 
lightcurves, the spectral properties should in principle be 
governed by the dominant pulse at a given time. The tail and 
high-latitude emission of a pulse can be observed only if it is 
followed by a period of inactivity. 

Examples A and B are very encouraging as they illustrate 
that low-energy slopes a can be expected in the range -3/2 < 
a < - 1 in the internal shock model with dominant synchrotron 
radiation. On the other hand, keeping the same assumption for 
the dynamics of the internal shocks as in case A and B, it seems 
difficult to find microphysics parameters leading to even steeper 
slopes — 1 < a < -2/3. This can be understood from the two 
shell model presented in the previous subsection. To reach the 
necessary condition r c /r m ~ 0.1 - 1, it is necessary to have 
collisions at lower radii, and/or with larger bulk Lorentz fac- 
tor, and/or to reduce the magnetic field (see Eq.fT2l>. This can be 
achieved in different ways : decreasing the contrast k, increas- 
ing the variability timescale t, increasing the Lorentz factor f , 
or reducing the kinetic energy flux E. In the following example, 
both the Lorentz factor and the kinetic energy flux have been 
changed : 

- Case C : the dynamics is the same as in case A and B ex- 
cept that the Lorentz factor has been multiplied by 3 and the 
kinetic energy flux reduced to E — 5 x 10 52 erg/s. The mi- 
crophysics parameters are =0.1 and ( — 10 -3 . This case 
is plotted in Fig. [10] and shows low-energy slopes a steeper 
than -1 during the whole duration of the pulse. The peak 
energy is £p i0 bs ~ 170 keV at the peak of the pulse. The ra- 
diative efficiency is still reasonably high (~ 60 %) however 
the end of the evolution occurs in slow cooling regime which 
results in a more complex evolution of the peak energy in the 
tail of the pulse than in the two previous examples. 

This example illustrates that the "marginally fast cooling 
regime" does provide low-energy slopes a > -1. However, fol- 
lowing Fig. [6] the conditions require a smaller radius and/or a 
low magnetic field. This tends to favor less energetic internal 
shocks. 

4.4. High-energy emission 

Interestingly, as already pointed out in IBosniak et al.l J2009), 
the scenario presented in this section - internal shocks with 
dominant synchrotron radiation in the soft gamma-ray range - 
require high Lorentz factors for electrons, which, because of 
Klein-Nishina corrections, always limits the efficiency of inverse 
Compton scatterings. So the high-energy spectrum (ferm//LAT 
range) does not show any bright additional component simul- 
taneously with the peak of the pulse in the GBM range, which 
seems in agreement with the GRB de tection rate of Ferm/ /LAT. 
However, as described in details in IBosniak et al.l d.2009), the 
physical conditions in the shocked region evolve during the in- 
ternal shock propagation in such a way that the w m parameter 
decreases in the tail of the pulse (see top panel in Figs. IMTOli. 
Inverse Compton scattering progressively enters the Thomson 
regime and becomes more efficient, especially in the low 
(high Fxh) case favored for steep low-energy slopes (see case 




Fig. 9. Case B : an example of a pulse generated by an inter- 
nal shock with synchrotron radiation in fast cooling regime 
affected by non negligible inverse Compton scatterings in 
Klein-Nishina regime. Same as in Fig. [8] The slope a is steeper 
(-3/2 < a < — 1 for the whole duration) than in case A, due 
to comparable values of w m but higher values of Yjh (see the 
theoretical interpretation in Sect.|3]l. 

B in Fig . [9j> . This leads to the delayed emergence of an addi- 
tional component in the high-energy spectrum. We will investi- 
gate in the future if this effect could explain the behaviour ob- 
served in ferm//LAT GRB lightcurves where delays between the 
GeV and the keV-Me V emission are ind eed observed (see for in- 
stance GRB 080916C. lAbdo et al.ll2009l) . 

5. Discussion and conclusions 

We present here a detailed discussion of the expected value 
for the low-energy slope a of the soft gamma-ray component 
(BATSE - Fenni/GBM range) in prompt GRBs, in the theo- 
retical framework where the kinetic energy of the outflow is 
extracted by internal shocks, and eventually injected in shock- 
accelerated electrons that radiate dominantly by synchrotron 
radiation. Our appr oach is based on the model developed in 
Bos niak et al.l(l2009l) . which takes into account both a full treate- 
ment of the dynamics of internal shocks and a detailed radiative 
calculation. 

1 . We show that in a large region of the parameters space of the 
internal shock model, the physical conditions in the emit- 
ting regions allow a combination of synchrotron radiation 
peaking in the soft gamma-ray range, and moderately effi- 
cient inverse Compton scatterings in Klein-Nishina regime. 
Interestingly, these scatterings affect the spectral shape of the 
synchrotron component, due to a better efficiency for low 
frequency photons. This results in a s teepening of the low- 
energy photon index a, with a — > - 1 (iDerishev et al.ll200Tt 
Bosniak et al. 2009; Nakar et al. 2009). For a large range of 
parameters regarding the dynamics of internal shocks (vari- 
ability timescale, bulk Lorentz factor, amplitude of fluctu- 
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a bright synchrotron component in the soft gamma-ray range, 
with a steep slope low-energy photon-index -3/2 < a < -1, 
as observed in most prompt GRB spectra. In this scenario, 
the additional Inverse Compton component at high energy is 
only sub-dominant and its intensity is not correlated to the 
intensity of the soft gamma-ray component. This seems in better 
agreement with FermifLPR observations and GRB detection 
rate than other scenario, such as standard SSC in Thomson 
regime, where bright components are predicted at high energy. 

On the other hand, even if steeper slopes in the range 
-1 < a < -2/3 can be achieved in the "fast marginally 
fast cooling regime", the corresponding constraints on the 
parameters of internal shocks are very strict. This regime could 
work if these steep slopes are associated, on average, with 
weaker intensities in the light curve. Indeed, the marginally fast 
cooling regime leads not only to a moderate radiative efficiency 
/rad — 0.5 - 0.7 (compared to / ra d — > 1 in the standard case), but 
also tends to occur in less energetic collisions. This limitation 
implies that this regime is probably not a robust explanation for 
all spectra with such steep slopes. 



Fig. 10. Case C : an example of a pulse generated by an 
internal shock with synchrotron radiation in "marginally 
fast cooling". Same as in Fig. [8] and [9] The slope a is steeper 
(-1 < a < -2/3 for the whole duration) than in case A and B, 
due to a higher ratio F c /F m . This ratio is of the order of 0. 1 - 1 at 
the peak of the pulse ("marginally fast cooling"), and even larger 
in the tail (slow cooling). 

ations of the Lorentz factor, kinetic energy flux), we pro- 
duce synthetic pulses with peak energies of a few hundred 
keV in the observer frame, and steep slopes in the range 
-3/2 < a < -1, at the peak of the lightcurve. The examples 
presented in the paper not only show high peak energies and 
steep slopes at the peak of the lightcurve, but also a general 
hard-to-soft spectral evolution, in agreement with observa- 
tions. This scenario constrains the microphysics in mildly 
relativistic shocks (shock acceleration and magnetic field 
amplification) : a large (e e ~ 0.1-1/3) fraction of the dis- 
sipated energy should be preferentially injected into a small 
(£ S= 0.01) fraction of electrons to produce a non-thermal 
population with high Lorentz factors, and the fraction of the 
energy which is injected in the magnetic field should remain 
low (ee ^ 10~ 3 ) to favor inverse Compton scatterings. The 
current knowledge of the microphysics in mildly relativistic 
shocks is unfortunately still rather poor and does not allow 
a comparison of these constraints with some theoretical ex- 
pectations. The prediction that only a small fraction £ of elec- 
trons are injected into a non-thermal power-law distribution 
leads naturally to a new question that we plan to investigate 
in the future : does the remaining quasi-thermal population 
of electrons contribute to the emission ? 
2. We also identify a regime of marginally fast-cooling syn- 
chrotron radiation with F c S5 F m which leads to even steeper 
slopes a — > -2/3 without reducing significantly the radiative 
efficiency of the process (/ ra d ^ 0.5). We present an example 
of a synthetic pulse in this regime, with a slope a > -1 for 
its whole duration. This requires low radii, and/or large bulk 
Lorentz factors, and/or weak magnetic fields. 

The present study shows that for a large region of the 
parameter space, internal shocks lead to spectra dominated by 



The scenario proposed in this paper can never produce 
low-energy slopes a steeper than -2/3, whereas such slopes 
are measure d in a non negligible fraction of GRB spectra 
dRvdel l2004t iGhirlanda et all l200l . About 20 % of GRBs 
have more than 50 % of their spectra with such very steep 
slopes. One possibility in such cases would be the appear- 
ance of a quasi-therma l component of photospheric origin 
dMeszaros & Reesl 120001: IDaigne & MochkovitchTl2002t iPe'erl 
2008; Belobo rodovl 120101: iPe'eretalJ teOloIT The emission 
from both the photosphere and internal shocks have a similar 
duration, the latter having only a very short lag behind the first. 
The intensity of the photospheric emission depends strongly on 
the unknown mechanism responsible for the acceleration of the 
relativistic outflow. A pure fireball would lead to a dominant 
thermal emission, whereas me chanisms suc h as the "magnetic 
rocket" recently suggested by iGranot et alj ((2010b) would on 
the other hand correspond to much colder jets with only a faint 
photospheric emission. In addition, the way it will superimpose 
on the non-thermal emission from internal shocks depends on 
the details o f the initial distribution of th e Lorentz factor in 
the flow (see IDaigne & Mochkovitchll2002l) . Very steep slopes 
could be observed in the time bins where the photospheric 
emission is dominant. 

Other effects could be a source of additional complexity in 
the prompt GRB spectrum within the scenario proposed in this 
paper. As it is required that only a fraction of available electrons 
are shock accelerated into a non-thermal power-law distribution, 
possible extra components in the spectrum could be associated 
with the remaining thermal population of electrons. Preliminary 
investigations show that the synchrotron radiation f rom these 
electro ns is easily self-absorbed, in agreement with IZou et all 
(120091) . On the other hand, due to lower electron Lorentz fac- 
tors, inverse Compton scatterings by these electrons usually oc- 
cur in Thomson regime, which is efficient. This could lead to 
additional components at low and/or high energy that offer in- 
teresting perspectives for the interpretation of the complex be- 
haviour observed in GRBs detected by Fermi[LPC[. We also find 
that in some cases, a non negligible fraction of the radiated en- 
ergy is re-injected in pairs due to yy annihilation. These pairs 
can radiate and scatter photons as well. These second order ef- 
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fects are not included in the present version of the model and 
could impact the final spectra shape. 
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